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Negative Magnetoresistance in the Superconducting State 
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We report on the c-axis resistivity p c (H) in Bi2Sr2CaCu2 0s+ < 5 that peaks in quasi-static magnetic 
fields up to 60 T. By suppressing the Josephson part of the two-channel (Cooper pair/quasiparticle) 
conductivity a c (H), we find that the negative slope of p c (H) above the peak is due to quasiparticle 
tunneling conductivity cr q (H) across the C11O2 layers below H C 2. At high fields (a) cr q (H) grows 
linearly with if, and (b) p c (T) tends to saturate (a c / 0) as T — > 0, consistent with the scattering 
at the nodes of the d-gap. A super linear cr q (H) marks the normal state above T c . 
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A clue to the mystery of high temperature super- 
conductivity is likely to arrive from the most peculiar 
normal-state properties of cuprate superconductors [Q. 
The observed 'pseudogap' features in the quasipar- 
ticle (QP) spectrum of the underdoped cuprates, such 
as Bi 2 Sr 2 CaCu208+5 (Bi-2212) |,§ and YBa 2 Cu30 7 -5 
[|],||, are taken by many as a signature of a non- Fermi- 
liquid electronic structure above T c 0. On the other 
hand, recent c-axis tunneling [|| data in the supercon- 
ducting state of Bi-2212 suggest that the <i-wave Fermi- 
liquid character of QPs is restored at low temperatures 
and low magnetic fields. The question arises where a 
crossover to the non-Fermi liquid behavior will occur and 
what is its signature in the c-axis transport - a direct 
probe of QP tunneling. Energies of the order of the gap 
Ao can be accessed with temperature, applied voltage, or 
magnetic field. The effect of field may be unique, since 
it can affect other degrees of freedom (e.g., magnetic ex- 
citations) that are coupled to QPs. 

The field dependence of the c-axis resistivity was first 
considered by Briceno et al. and by Gray and Kim || 
in the context of a 'giant' magnetoresistance vs tempera- 
ture peak observed in Bi-2212. They ascribed the p c (T) 
peak to two competing processes: (1) Josephson and (2) 
quasiparticle tunneling, the latter masked by supercon- 
ductivity below T c . The former, understood only recently 
by Koshelev Jl0|| , leads to a positive c-axis magnetoresis- 
tance at sufficiently low fields |9 |Tl|JT^ |. At high fields, the 
observed crossover to a negative slope in p c (H) has been 
attributed to the normal state [ |TT|— p^4[ and to the emer- 
gence of the pseudo- or spin-gap |12|]. Consequently, the 
field at the peak in p c (H) in Bi-2212 was taken as L H C 2 
[|l4f , in an apparent conflict with the if-dependencies of 
p ab (H) and in-plane T c2 (H) §,0Q. 

In this Letter we resolve the origin of negative c-axis 
magnetoresistance in layered high-T c superconductors by 
exploring quasiparticle dissipation in magnetic fields up 
to 60 Tesla. We demonstrate that at high fields, above 



the maximum in p c (H) at iJ*, negative magnetoresis- 
tance in Bi-2212 is controlled by QP tunneling in the su- 
perconducting state. We access QP current by suppress- 
ing the Josephson current in two ways: (i) with transport 
current in thin Bi-2212 mesas, and (ii) with a 60 T field, 
which in underdoped Bi-2212 crystals exposes QP tun- 
neling down to ~ 22 K. We show that high- field c-axis 
conductivity is linear-in-field, in agreement with the field- 
linear QP conductivity <j q we find in Bi-2212 mesas in the 
resistive state. Near 60 T, at low temperatures p c tends 
to a finite (saturation) value, as predicted for a d-wave 
superconductor. The crossover to the normal state above 
T c is witnessed by a superlinear a q (H), plausibly related 
to the pseudogap. 

We used high quality, optimally and slightly under- 
doped Bi-2212 crystals (T c 's « 92.5 and 89 K respec- 
tively). Here we will show data on the underdoped crystal 
(900 x 600 x 20 /im 3 ), since it afforded us a downward- 
expanded temperature range of negative magnetoresis- 
tance. The resistivity was measured in magnetic fields 
along the c-axis up to 60 T generated in the Long Pulse 
(LP) System at the National High Magnetic Field Labo- 
ratory (NHMFL) in Los Alamos, NM. One advantage of 
the LP system is that it allows us to avoid sample heat- 
ing by the induced eddy currents - the main problem 
encountered in short-pulse experiments, where a large 
dB/dt > 10 4 T/s produces significant heating and associ- 
ated thermal hysteresis effects for larger than micron-size 
samples. The LP system delivers a nearly triangular field 
pulse with the fastest up- and down-ramp rates ~ 150 
and 360 T/s and a duration ~ 2 sec. A non-hysteretic 
resistivity at these different rates assures a minimal ef- 
fect of eddy currents. The temperature was controlled to 
better than 50 mK, with the massive copper-dust /epoxy 
sample holder serving as a thermal anchor. A standard 
4-probe contact configuration was used to measure resis- 
tance by a lock-in technique at 17 kHz, recording it with a 
fast analog-to-digital converter. The current- voltage (I- 
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V) characteristics of mesa-patterned high quality Bi-2212 
whiskers, prepared by a Focused Ion Beam technique [||, 
were measured at NHMFL in Tallahassee, FL. 

Figure 1 shows the raw data of the c-axis resistivity p c 
vs magnetic field for the Bi-2212 crystal in the tempera- 
ture range 22.5K - 110K. We point to two major features 
below T c = 89 K: first, a maximum in p c (H) and second, 
the near saturation of high-field p c vs T at the lowest T 



(inset in Fig. 1) |15|. The latter is a new observation. 
Note that above 35 K, p c (T) at 55 T roughly follows a 
InT dependence up to ~ T c . Above T c , p c has a weaker 
(~ 0.1%/T at 110 K) field dependence. Remarkably, a 
60 T field even at 70 K does not suppress p c to it's value 
above T c . 
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FIG. 1. Out-of-plane resistivity p c vs H of Bi-2212 crystal 
in fields up to 60 T at different temperatures, j = 0.05 A/cm 2 . 
Inset: p c vs T at 55 T (♦) and at T (line). 

A maximum in p c (H) at H*(T) is summarized in the 
if — T diagram in Fig. 2. Below if* the magnetoresis- 
tance is positive and above negative. Such independence 
has been observed before []T2||L4| , and attributed to a 
crossover to the normal state. Indeed, since there is 
a (weak) negative magnetoresistance above T c (Fig. 1), 
one may easily fall into the trap of identifying if C 2 with 
if*. However, p c is a measure of tunneling across the 
Cu02 layers and a direct signature of the normal state 
should come from the in-plane (non-tunneling) resistiv- 
ity p a b(H). It reveals a maximum g|| at if > if* (inset 
in Fig. 2): at if*, p a b has still a strong positive magne- 
toresistance, clearly originating from the superconduct- 
ing state. So, this is a first telltale sign that if C 2 > if*, 
since the changes in p a b(H) and p c (H) do not track. Thus 
we take the maximum in p a b(H) as an underestimate of 
H c2 H (see Fig 2), the assignement we will confirm in 
the considerations that follow. 

Below we show that a key to understanding magnetore- 
sistance is the realization that the c-axis conductivity is a 
parallel, two-channel tunneling process: (i) aj of Cooper 



pairs (mostly at low fields) and (ii) a q of quasiparticles 
(dominant at higher fields): <j c = aj + <j q . We first 
discuss a j. Consider a stack of superconducting sheets 
spaced by s subject to a c-axis magnetic field. Dissipa- 
tion for Josephson interlayer current is caused by phase 
difference slips due to motion of the Josephson strings 
attached to pancake vortices in adjacent layers (To). Dif- 
fusive drift of pancakes governed by pinning and thermal 
fluctuations results in motion of current-driven Joseph- 
son strings. The barrier for pancake motion is large at 
low fields but lessens with increasing; field §. The field 
dependence is explicit in the derived universal relation- 
ship between the in-plane and the c-axis conductivity M 



aj a ab ^ s 2 Ej/(BT 2 ) (x a ab /B, 



(1) 



where Ej is the Josephson energy per unit area at B = 0, 
B ~ H, and cr ab ex exp(/7 /T) is controlled by the energy 
barriers U(H) to thermally activated pancake hopping 
0. For a 2D vortex lattice U(H) ~ -lnif flf], implying 
a power-law field dependence of aj = aH~ v with v(T) w 
1 + U/T. The low-field (J a b/&c is indeed nearly linear in 
if (inset in Fig. 3), confirming Eq. 1. 
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FIG. 2. H - T diagram showing H*(T), the location 
of peak in p c (H) (•). Rough estimates of H C 2(T) from the 
broad maximum in p a b(H) |l2| (A), from the demise of the 
j-dependence of p c (H) (♦, see Fig. 3), and from 1/(3 (▼) in 
Eq. 2. Inset: p c and p a b vs H. A positive dp a b/dH persists 
well above iT(T). 

The Josephson nature of <jj implies that this contri- 
bution and thus if* may be significantly affected by an 
injection of extra interlayer transport current j. This 
current - in addition to pancake disorder - suppresses 
Josephson coupling, while its effect on a q is negligible. 
For the data of Fig. 1, a downshift of if* at larger j is 
evident in Fig. 3. At a higher temperature (55 K) the j- 
dependence of p c (H) (i.e., Josephson current) smoothly 
disappears above 45 T. This crossover into an ohmic dis- 
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sipation regime - controlled mainly by the QPs []l9| - 
further marks a lower bound to H C 2 > if* (Fig. 2). 

An immediate consequence of this picture is that at 
high fields aj <C cr g , and thus a c ~ a q . To unambiguo- 
suly separate the field behavior of the two channels, we 
first independently test the field dependence of <j q via c- 
axis I-V characteristics in a tiny (~ 2/im 2 ) step-like mesa 
bar (sketch in Fig. 4) carved out of a Bi-2212 whisker, 
with the thickness of ~ 50 Cu02 layers. Here, the Joseph- 
son current can be suppressed by a fairly low transport 
current || and the decreasing branch of the I-V curves 
in the resistive state is determined purely by quasiparti- 
cle tunneling. The inset in Fig. 4 shows a set of f-Vs 
from which QP conductivity (the initial slope of I vs V) 
was extracted as in Ref. ||. The obtained a q (H,T) vs 
c-axis DC field is remarkably linear up to 33 T at all 
temperatures (main panel of Fig. 4). 
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We turn now to the entire field and temperature depen- 
dence of the c-axis conductivity in a macroscopic crystal 

* C (H, T) * oJT^ + (r g (0,0)fry + /?H] ; (2) 

Fig. 5 shows <j c (H) for the crystal. Below T c there is 
clearly a decreasing aj below if* and, for if ^> if*, the 
if -linear cr q . Both channels are well described by Eq. 2 
as illustrated by a fit at 55 K. Such fits give v in the 
range 1.5 - 3.5 between 70 and 22.5 K, confirming the 
In if dependence of hopping barriers with U ~ 45 K at 
low T. The coefficient a(T) increases with decreasing 
T [|l7j. The temperature dependence of the zero-field 
<jg(0, T) gives rj(T) = 1 + cT 2 (inset in Fig. 5), in agree- 
ment with the results on mesas. From the coefficient 
c = 7r 2 /187 2 (T < 7) H we find the effective scattering 
rate of QPs due to impurities inside layers, 7 ~ 0.4T C . 



And, significantly, we obtain a nonzero extrapolated 
a q (H -> 0,T 0) [« 2.5 (kfi cm) -1 ]. This result is a 
signature of a d-wave superconductor Q an d it confirms 
direct measurements of cr g (0, 0) ~ 1.5 — 3.7 (kQ cm) -1 in 
mesas Q . This explains the saturation of high-field p c at 
T — > (Fig. 1) which naturally reflects the nonzero <j g 
arising from the scattering by impurities in the nodal re- 
gions of the gap. Above T c , the two-channel description 
of Eq. 2 obviously fails. The o~ q (H) becomes weaker and 
is no longer linear (Fig. 5). The change across T c is very 
gradual and subtle, consistent with the existence of the 
pseudogap above T c [|],[|. 




10 20 30 

1^ h CD 

FIG. 4. Normalized quasiparticle c-axis conductivity as a 
function of H || c obtained from the I — V curves (top inset) 
measured on the mesa-shaped Bi-2212 (sketched). 

We now address our result - the linear increase of <j q 
with H. Within the Fermi-liquid model, cr q (H) in the 
superconducting state is proportional to the QP density 
of states (DOS) and inversely proportional to the effec- 
tive scattering rate for interlayer tunneling, j c . In the 
s-wave picture, a q w ctnH/H C 2 (<tat is the normal state 
conductivity), since the DOS is proportional to the num- 
ber of vortex cores. However, in a d-wave superconductor 
the dominant contribution to a q comes from the fourfold- 
symmetry nodal protrusions of delocalized states outside 
the cores |2l[ |. They also lead to the increase of DOS 
with H due to the Doppler-shifted QP energy caused by 
the in-plane super currents around vortices. Recent cal- 
culation Q shows that to leading order in the Doppler 
shift energy with respect to Ao, the increase of DOS is 
compensated by an increase of j c . The increase of j c is 
caused by random positions of pancake vortices result- 
ing in different Doppler shifts between equivalent points 
in adjacent layers. This compensation results in a field 
independent cr q (H,T) w ^(0, 0)r?(T) in Eq. 2. But, at 
high fields the core contributions and other corrections 
(i.e., nonlinear QP spectrum near gap nodes |2^| ) must 
enter. These corrections give a net increase of cr q (H) 
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as - a q (0,0)H/H A . H A = ^/Sl/h 2 v 2 F is « 40-80 T 
at low T |2^| . In the BCS theory if a corresponds to 
H C 2 ~ $0/(^0 2 and 1/^8 (Eq. 2) gives a rough estimate of 
H C 2(T) plotted in Fig. 2. However, this estimate is fuzzy 
since the above does not include the pseudogap which 
may complicate the magnetoresistance at high fields. 

Note that at a 60 T field a q (H, T) is only w 2a q (0, T), 
far below the normal state value: e.g., at 140 K <j c ~ 
55(k^ cm) -1 ~ 20<j g (0,0). It is also below a at reached 
by applying voltage V > 2Ao/e at low T p|,p3[. This is 
surprising. It may reflect the influence of the pseudogap 
which has been shown to be robust to changes in the 
magnetic field ||. This is consistent with the gapped 
QP spectrum inside vortex cores observed by STM H. 



tract No. AL99424-A009 and by LANL LDRD Grant 
No. 98037. The work of Yu.I.L. was supported by Japan 
Science & Technology Corp. and by the Russian State 
Program on HTS, Grant No. 99016. 




FIG. 5. a c (H) below and above T c = 89 K. A fit 
at 55 K to a superposition of Cooper pair (dash) and 
quasiparticle (dash-dot) contributions in Eq. 2 is indicated. 
The high- field a c (H) becomes nonlinear above T c . Inset: 
Zero-field a q vs T/T c extracted from a c (H,T) for a Bi-2212 
crystal for j = 0.05 A/cm 2 (•) and j = 0.1 A/cm 2 (y) 5 
and obtained from the I-V curves for a mesa (o). Both fit a 
independence (dash) up to T ~ 7. 

In summary, the field dependence of the c-axis conduc- 
tivity in Bi-2212 is consistently described by tunneling of 
Cooper pairs and of quasiparticles across a d-wave su- 
perconductor. We demonstrate that the maximum in 
p c (H) comes from the competition between these two 
conduction channels below H C 2 . The low-T saturation of 
the high-field p c (T) is a consequence of the d-wave char- 
acter of the superconducting state. Remarkably, cr q (H) 
remains linear in H and small up to 60 T. It becomes 
gradually weaker and superlinear above T c , where the 
QP tunneling is controlled by the pseudogap. 
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